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It is shown that the critical velocity of a small sphere oscillating in superfluid helium increases with 
the square root of the oscillation frequency. This behavior can be described by a simple dimensional 
argument. The size of the sphere and the temperature of the superfluid are found to have no or only 
very little effect. Surface properties of the sphere and remanent vorticity may have an influence but 
have not been under systematic investigation in these measurements. 
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A quantitative description of the critical velocity for 
the onset of turbulent drag of a macroscopic body moving 
in superfluid helium is still an unsolved problem although 
there is a large amount of experimental data available 
from vibrating wires, grids, quartz tuning forks, and os- 
cillating spheres. At small driving forces F the velocity 
amplitude v grows linearly because the drag force due 
to quasiparticle scattering is proportional to the velocity. 
At a certain critical value v c vorticity is shed by the mov- 
ing body and the drag increases strongly. Typically, these 
velocity amplitudes are in the range from mm/s up to al- 
most m/s and therefore much lower than the Landau 
critical velocity of about 50 m/s. Above v c the velocity 
amplitude grows approximately with the square root of 
the driving force because turbulent drag is proportional 
to the square of the velocity Q. Critical velocities for 
vortex production have been calculated, for a review see 
3]. For example, a vortex ring of radius R is calculated 
to expand if a critical velocity is exceeded that is given 
approximately by Q 
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where n w 10~ 7 m 2 /s is the circulation quantum and 
ao ~ 0.15 nm is the size of the vortex core or the co- 
herence length. For a sphere of radius R = 100 /im this 
would give 2.5 mm/s for a ring of the same size. This 
is at least an order of magnitude lower than the experi- 
mental value, see below. Within this picture, it follows 
that smaller vortex rings will be shed, probably because 
of the roughness of the surface. Numerical simulations 
of vortex shedding have been performed using full Biot- 
Savart calculations and assuming a smooth sphere [l], [B| . 
Numerically, the obtained critical velocities are of same 
order as in experiments. However, the results are some- 
what sensitive to the computational accuracy and there- 
fore the frequency dependence for v c is not extracted yet. 
In any case, a simple dimensional argument tells us that 
a characteristic velocity scale in the turbulent superfluid 



must be given by k/R, whatever the length scale R is. 

The question arises whether the physics will be dif- 
ferent in case of oscillatory motion. The turbulent drag 
force Fd — "/v 2 ~ Fo (7 = CdP^R 2 /2, where p is now the 
superfluid density and C'd ~ 0.4 for a sphere @) differs 
from the classical result for steady flow only by the shift 
Fo • This shift implies a critical velocity for the on- 
set of turbulent drag given by v c — t/Fq/j. The shift is 
temperature independent Q , and so is this critical veloc- 
ity. A theoretical calculation of Fo is not available. (It 
should be mentioned, however, that a similar drag force 
has been calculated for a body moving through a Bose 
condensate Q.) In this situation it is necessary to resort 
to dimensional considerations. Because there is a time 
scale now given by the oscillation frequency u>, we have a 
new velocity scale for the onset of the turbulent regime 
given by Q 



(2) 



Remarkably, there is no length scale present any more. 
There might be, however, a logarithmic term like the one 
in Eq.(l). This implies that this critical velocity docs 
not or only very little depend on the size of the body. 
In the following, this scaling law is confirmed experimen- 
tally with spheres oscillating at various frequencies. In 
Appendix 1 a model calculation based on Kelvin waves 
is presented that supports the scaling law. In the Ap- 
pendix 2, the same dimensional arguments are applied 
to a classical viscous liquid and shown to be consistent 
with analytical calculations for spheres. 

The experiments were performed at Regensburg Uni- 
versity in the years from 1993 to 2000. Most of the results 
have been published in earlier work (Toj , but a systematic 
analysis of the frequency dependence of the critical veloc- 
ities was only performed recently. The data were taken 
with 2 different spheres, having a radius of 100 /im and 
124 /im. The smaller one was used in a 3 He cryostat, the 
other one in a dilution cryostat. The magnetic spheres 
were placed inside a horizontal parallel plate capacitor 
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consisting of niobium electrodes to which a voltage was 
applied while the temperature was dropped below the su- 
perconducting transition temperature of niobium, ca.10 
K. At helium temperatures the charged spheres were lev- 
itating between the capacitor plates. Frozen flux in the 
electrodes was necessary for a stable vertical oscillatory 
motion when a small ac voltage was applied at resonance. 
The frequency of the oscillations could not be controlled 
quantitatively because it depended on the frozen flux of 
the particular levitation status. It was observed, how- 
ever, that a faster cooling rate from 10 K to 4.2 K in 
general lead to higher resonance frequencies. The criti- 
cal velocities were measured by increasing the driving ac 
voltage from the linear flow regime to the turbulent one 
and back. The results of the various runs, all taken below 
1 K, for both spheres are displayed in the Figure 1. 
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FIG. 1: Critical velocity for the onset of turbulence as a func- 
tion of the frequency of 2 oscillating spheres. Blue squares: 
radius 124/im; red dots: radius 100/im; green solid line: 
v c = 2.85 V K2nf (m/s); temperature is below 1 K. 



thermally activated process for the onset of turbulence 
in these experiments. 

It will be interesting to compare these results for 
spheres with those for vibrating wires. Earlier observa- 
tions by the Osaka group indicated frequency indepen- 
dent critical velocities below 2 kHz (Tj|. A more recent 
analysis by other authors did show an increase scal- 
ing approximately with w 1 / 3 . A more complete and very 
recent work of the Osaka group is in fact in agreement 
with Eq.© Q. 
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APPENDIX 1 

One simple way to derive the scaling relation for v c 
(Eq.©) is the following. The dispersion relation for 
Kelvin waves with wave number k is given by: 

If we assume that the external drive at frequecy to only 
drives the scales R = 2ir/k given by the above dispersion 
relation, we find 



It is obvious that the predicted y/nuj variation gives a 
fair fit to the data. The errors of the individual values 
of v c vary from ca. 1% to about 5%, depending on how 
closely v c was approached in the experiments. The large 
scatter visible near 200 Hz is by far larger than any error 
bar. Recently, it has been found that remanent vortic- 
ity affects the critical velocity of vibrating wires and its 
absence greatly increased v c [Tlj] . This might be a rea- 
son for the scatter. Moreover, surface properties might 
affect a numerical prefactor in Eq.©, although, as indi- 
cated above, only on some logarithmic scale, otherwise we 
would have had much more scatter. The mostly higher 
resonance frequencies of the smaller sphere are due to its 
smaller mass and probably also due the faster cooling rate 
from 10 K used with the 3 He cryostat. Towards higher 
temperatures only a rather weak increase of v c by about 
10% was observed up to 2 K [l2] |. This rules out any 



R w ^ktt ln(2/fca )/w. (4) 
Inserting this then in Eq. (pj gives a critical velocity 

ln(167r/fca ) , — ,_. 

v c s» =y/-ca;. (5) 

2tt^/it ln(2/fca ) 

The coefficient in front should in principle be calculated 
self consistently using the dispersion relation but depends 
only very weakly on the frequency. For typical sphere ex- 
periments (200 Hz) the coefficient should have a value of 
0.40 which is somewhat smaller than 2.85 obtained from 
the data of Fig.l. Of course, there is no real proof that 
Eq. ([!} should give the correct value in case of oscillating 
flows. 
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APPENDIX 2 



The same dimensional arguments are applied now to 
a classical viscous liquid. For the stationary case the 
velocity scale is set by v/R, where v is the kinematic 
viscosity. This is confirmed by comparing the laminar 
drag force of a sphere Aw with the turbulent drag force 
jv 2 . From Stokes' solution we have A = QirvpR. Hence, 
by extrapolating both regimes, the drag forces are equal 
at a characteristic velocity 
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For the oscillating case we get from dimensional argu- 
ments a velocity scale y/vui. Stokes' solution for the lin- 
ear drag is now 6] A = Q-rtvpR- (R/5), where S = \j2vjuj 
is the viscous penetration depth, and it is assumed that 
R 5, which usually is the case for macroscopic bod- 
ies in liquid helium at or above the superfluid transition 
temperature. The velocity scale is now given by 
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Note that there is no length scale in this expression and, 
therefore, v c h does not depend on the size of the spheres. 

Comparing these results with those for the superfluid, 
we note the similarity, except that the kinematic viscos- 
ity v is replaced by the circulation quantum k. The ma- 
jor differences, however, are, firstly, that in the classical 
liquid there are about three orders of magnitude in ve- 
locity between the Stokes regime of linear drag and fully 
developed turbulence and, secondly, because of the ex- 
trapolation, v c h is not a point on the v(F) response of 
the body, whereas in the superfluid we have a true crit- 
ical velocity of the oscillator for the onset of turbulent 
drag. Experiments with quartz tuning forks, for which 
analytical expressions of the drag forces are not available, 
confirm this classical scaling 
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